Entanglement and level crossings in frustrated ferromagnetic rings 
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We study the entanglement content of a class of mesoscopic tunable magnetic systems. The 
systems are closed finite spin-1/2 chains with ferromagnetic nearest-neighbor interactions frustrated 
by antiferromagnetic next-nearest-neighbor interactions. The finite chains display a series of level 
crossings reflecting the incommensurate physics of the corresponding infinite-size chain. We present 
dramatic entanglement signatures characterizing these unusual level crossings. We focus on multi- 
spin and global measures of entanglement rather than only one-spin or two-spin entanglements. 
We compare and contrast the information obtained from these measures to that obtained from 
traditional condensed matter measures such as correlation functions. 



I. INTRODUCTION 

The study of entanglement measures in the eigenstates 
of condensed matter Hamiltonians is a rapidly growing 
field of interdisciplinary research [ij . While a great many 
results have been reported, the typical situation is that 
entanglement studies of condensed matter systems pro- 
vide alternative signatures of phenomena already known 
from more traditional condensed-matter techniques. In 
this article, we use measures from quantum information 
to characterize a spin system which is poorly under- 
stood in the condensed matter literature, namely, spin 
rings with ferromagnetic nearest-neighbor interactions 
frustrated by antiferromagnetic next-nearest-neighbor in- 
teractions. We show how the unusual properties of this 
system are reflected in quantum information quantities. 

The systems we present are closed rings of 0(10) 
spins. The study of finite-size quantum systems have en- 
joyed a huge resurgence in the past decade because of 
cold-atom and nanoscience developments, which provide 
the context for designing and studying finite-size Hamil- 
tonians for their own fundamental properties rather than 
as mere tools for investigating the thermodynamic limit 
N ^ OO. We will focus on a sequence of ground-state 
level crossings in such systems, which reflect lattice- 
incommensurate correlations. We show that certain mea- 
sures of entanglement display remarkably strong param- 
eter dependences near these level crossings. Our calcu- 
lations of the concurrence show that the entanglement 
structure of the states between successive level crossings 
is ideal for tuning the distance between spins which are 
mutually entangled. We also present some results on en- 
tanglement dynamics, again focusing on the effect of the 
level crossings. 

It is natural to think of ground state level crossings as 
finite-size analogs of first-order phase transitions in the 
thermodynamic limit. As a corollary, one expects that 
ground state properties change suddenly at the crossing 
('transition') point and that away from the exact crossing 
point the system possesses no knowledge of the proximity 
of the crossing. In this regard, our quantum information 
quantities reveal a peculiarity of the crossings in our spin 
rings, because they seem to evolve in response to the 



proximity of the level crossings, already at some distance 
away from the parameter value where the crossing oc- 
curs, i.e., the ground state is "aware" of the nearby cross- 
ing. This counter-intuitive phenomenon is clearly seen in 
multi-spin (but bipartite) entanglements as well as fideli- 
ties, but not in traditional condensed-matter quantities 
such as correlation functions. Note that, since the level 
crossings are not expected in the thermodynamic limit, 
this issue is meaningful only in our mesoscopic context. 

To calculate a bipartite entanglement, one has to first 
decide on a partition. For a condensed-matter model or 
even its finite-size versions, there are a large number of 
ways in which the system can be bi-partitioned. A variety 
of different partitionings, and correspondingly different 
entanglements, have been used for entanglement studies 
in condensed matter models. Many studies have focused 
on the entanglement between two sites as measured by 
the concurrence or negativity {e.g., [3, d, 0, H, 0, [1]), 
or the entanglement between one spin/site with the rest 
{e.g., 9, 10, 11]). Such measures generically are not ex- 
pected to provide information different from traditional 
condensed matter quantities like 2-point correlation func- 
tions, since the entanglement measure can be related eas- 
ily to correlation functions {e.g., [1, H, 0, S ) ■ In this 
sense, many-site measures are more innovative. Exam- 
ples are block entanglement in ID and in 2D [l^ , and 
more recently studied spatially disconnected partitions 
such as sublattice entanglement il4| and entanglement 
between labeled itinerant particles Such entangle- 

ment measures are not easily related to usual correlation 
functions, and therefore hold the possibility of providing 
genuinely new tools for probing collective phenomena in 
many-particle systems. For our finite-sized chains, we 
will present entanglement results for various partitions, 
both two-site and many-site. 

The dynamics of entanglement measures in many- 
particle systems has attracted some attention recently 
[l6f|, but is still a relatively poorly understood topic, re- 
fiecting the general difficulty of non-equilibrium physics 
in quantum condensed matter. Entanglement dynam- 
ics has also gained interest from the quantum informa- 
tion perspectivCje.g., in the context of adiabatic quan- 
tum computing [l7| and in the context of quantum state 
transfer ('quantum communication') through spin chains 
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FIG. 1: (a-c) Various partitions of closed chain for calculating 
bipartite entanglement entropy. Spins belonging to partition 
A (B) are represented by open (full) circles, (d) Geometry 
for exploiting the sublattice entanglement Ss.i., e.g., in an 
experimental realization. In this setup Ss.i. is simply the en- 
tanglement between inner and outer rings. Dashed and full 
lines represent Ji and J2 interactions respectively. 



|18| . Obviously, the number of different conceivable non- 
equilibrium situations is large; we are forced to restrict 
the types of dynamics considered. In this article we will 
describe the entanglement dynamics after a sharp change 
{quench) of interaction parameter. 

After describing in Sec. |TT] the various entanglement 
measures and partitions to be used, in Sec. IIIII we intro- 
duce the Hamiltonian and its condensed-matter context. 
Sec. IIIII also describes the level-crossings and compares 
to other examples of similar crossings known to us. The 
next sections present entanglement results, starting in 
Sec. IIVI with the quantities containing the most tradi- 
tional information, namely, the concurrences which are 
readily related to correlation functions. Sections |V] and 
IVII present quantum-information quantities which are 
more distinct from condensed-matter measures. Sec. IVIII 
describes our results on entanglement dynamics. 



II. PARTITIONS AND ENTANGLEMENTS 

We will mostly measure entanglement using the en- 
tanglement entropy. The entanglement entropy between 
partitions A and B is defined using the reduced density 
matrix of one part (e.g., pA — tr^ p obtained by tracing 
over B degrees of freedom) to calculate the von Neu- 
mann entropy, Sa = — tr[pAliipA]- Since the purpose of 
using entanglement in many-particle systems is to pro- 
vide viewpoints not present in usual condensed-matter 
measures, we want to consider more 'non-local' entan- 
glements. The partitions we choose for this purpose are 
shown in Figure [TJ 

The most obvious choice is for A and B to be spatially 
connected blocks, in particular for each partition to be a 
connected half of the ring (Fig.[T^). This partitioning is 
suggested by many existing results on block entanglement 
entropies in one dimension [12]. 

Another choice is to take all the odd sites as parti- 
tion A and the even sites as partition B — this is then 
the entanglement between A and B .sublattices, in usual 
quantum magnetism language (Fig. [T]d). The entangle- 
ment between sublattices has been studied for a number 



of spin chains and lattices, and extrema of Sa have been 
claimed to characterize some quantum phase transitions 
Note also that since it is natural to think of the 
J1-J2 chain as a 'zigzag ladder' {e.g., [H, [13, HH), the 
sublattice entanglement is simply the entanglement be- 
tween the two legs of such a ladder. In the mesoscopic 
(finite-size) context that we are interested in, a J1-J2 ring 
is likely to be implemented as shown in Fig. [TJi; the two 
'legs' are now two J2 rings (either inner or outer rings, 
or an upper and a lower layer), connected by zigzag Ji 
interactions. From an experimental quantum informa- 
tion processing perspective, the sublattice entanglement 
is natural as it is simply the entanglement between inner 
and outer (or upper and lower) rings. 

Since the antiferromagnetic interaction in our case is 
the next-nearest-neighbor interaction, we extend the idea 
of the sublattice partition {ABAB...) and also use the 
AABBAA... partition shown in Fig.[TJ:. 

We call the entanglement entropies corresponding to 
the above three partitions respectively 81/2, Ss.i., 82/4- 

Another nonlocal quantity motivated by quantum in- 
formation theory is the fidelity, the overlap between 
ground state wavefunctions at slightly different system 
parameters [22]. We describe fidelities in our spin rings 
in Section IVII 

We also present results for more 'local' measures of 
entanglement, namely, (1) entanglement entropy for 2- 
site partitions, and (2) the concurrence Cr measuring the 
entanglement between sites i and i-|-r in the environment 
provided by the rest of the chain [23| . The concurrence 
is defined as 



Cr — max 



— V A2 — V A3 — 



where are the eigenvalues in decreasing order of the 
matrix pA{o'y ® <^y)PAi'^v ® (^y)! and pA is the reduced 
density matrix of the two-site subsystem. Since we will 
consider singlet states, two-site reduced density matrices 
are strongly constrained, and the concurrence is a simple 
function of the spin-spin correlation functions. Two spins 
are entangled {Cr > 0) if the correlations between them 
are sufficiently antiferromagnetic, {SfS^_^_j.) < —1/12, 
and in that case Cr = —6{SQSr) — \ (derived in, e.g., 
0; see also Ref. 22 in [21]). The behavior of concur- 
rence in various spin models is briefly reviewed from the 
literature in Sec. IIVI in order to provide context for our 
concurrence results. 



III. FRUSTRATED FERROMAGNETIC RINGS 

We are concerned with closed chains of N localized 
spin-1/2 objects, i.e., a finite-size spin chain, interacting 
via the Hamiltonian 

H ^ Ji Si • Si+i -I- J2 Si • Si+2 , (1) 

i i 

where Ji < J2 > 0, and i is the site index obeying peri- 
odic boundary conditions appropriate to a ring. We de- 
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fine (3 = J2/IJ1I and (3 = f3~^. This spin chain is known 
for any even N to have a singlet ground state for < 4 
and a (A^+ l)-_fold degenerate ferromagnetic ground state 
manifold for /3 > 4 [24!]. At /? = 4 a singlet ground state 
is degenerate with the ferromagnetic manifold; these de- 
generate states are known exactly [l^ H^, HE 0, US] • 

We will confine ourselves to the singlet region, /3 € 
(0,4). The numerical results are presented for ring sizes 
divisible by 4, in order to avoid complications due to the 
parity of N/2. The rings for other even sizes have very 
similar behavior. 

Although far less studied compared to the Ji > 
case (Majumdar-Ghosh chain), the frustrated ferromag- 
net Hamiltonian has already been known for some time to 
possess some peculiar properties [l^, [13, [H, [2^ [13, HH . 
More recently, this model has attracted a fresh sur ge o f 
theoretical attention [H, [H, [H [H [H, [13, [li,lli|, 
partly due to the appearance of materials whose s pin 
physics is described by this model (see references in [20l . 
[32|l. Most studies concentrate on finite magnetic fields 
or finite temperatures and on properties in the thermody- 
namic limit 20, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 3§, 
or very near the transition point at /3 = 4 [27|, (S^, |39|, (4^ . 

Proliferation of ground state level crossings — An un- 
usual feature of the singlet region (3 G (0, 4) is the number 
of ground state level crossings, which proliferate as we in- 
crease the number of spins [2^. For N = 8, 12, 16, 20,... 
spins, there are respectively one, two, three, four,... level 
crossings in this region. A very likely interpretation is 
that the natural correlations in the chain for f3 £ (0, 4) 
are incommensurate with wavevector dependent on /3. 
This incommensurability is resolved in different ways by 
the different ground states. For larger rings, there are a 
greater number of possibilities to accommodate the in- 
commensurate correlations into the finite system, and 
hence the larger number of level crossings. The ground 
state lattice momentum alternates between and tt in 
the regions between successive level crossings. 

Very similar level crossings appear in other spin models 
known or thought to have spiral correlations. The exam- 
ples known to us are (1) the Majumdar-Ghosh chain for 
J2 > 0.5Ji dH; (2) the 'sawtooth' chain of Ref. 
(3) the frustrated ferrimagnetic chain of Ref. [4^. This 
phenomenon of proliferating level crossings in finite-size 
rings as a result of incommensurate correlations, how- 
ever, is not widely known or discussed in the condensed 
matter literature. 



IV. CORRELATION FUNCTIONS AND 
CONCURRENCES 
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FIG. 2: Correlation functions and related quantities for = 
16. (A) correlation function at /3 = 2. (B) The structure 
factor, S{q), at P — 2. (C) Position of the maximum of S{q), 
as a function of /3. (D) Concurrences between two spins at 
distance r — 2, 3, 8. (E) Entropy of entanglement between 
partition consisting of sites (i, i + r), and the remaining N — 
2 sites. Shown also for r — 2, 3, 8. In (C-E), nonlinear 
horizontal scale highlights the region near P — 4. 



tion of r = |z — j|. Between /? = and the first crossing, 
the signs of this quantity has the structure 



for r — 1,2,3,4,...., indicating that the state in this re- 
gion has a 'spin-density-wave'-like Ittiitt • • • ) structure. 
An example is shown in Fig.[2j^ for A^ = 16. The oscillat- 
ing magnitude /sign pattern also suggests 'nematic'-like 
correlations [33l . |44| . |45|; detailed distinctions are inap- 
propriate since this is not a macroscopic phase. 

The sign pattern of {SqS^.) changes at each level cross- 
ing; for A^ = 16 the patterns in the four regions (sepa- 
rated by level crossings) arc 



In this section we present more local quantities, in 
particular traditional condensed matter quantities like 
spin-spin correlation functions and structure factors, and 
closely related entanglements like concurrences. 

For singlet states, the correlation function (Si • Sj) is 
three times (S^Sj); we will discuss the latter as a func- 



The structure factor S{q), which is the Fourier trans- 
form of the real-space correlation function, has peaks at 
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<?max = {2Tr /N)m, with the integer m decreasing in unit 
steps at each level crossing, from m — N/2 at (3 = 
to m = 1 near the ferromagnetic transition point [2^ . 
Similar behavior (Fig. is also seen in the other ex- 
amples known to us where finite-size spin systems have 
level crossings due to spiral correlations 4l|, 43 |. 

We now consider the concurrences Cr between spins 
i and i -I- r on the rings. For most simple spin models, 
the concurrence does not extend significantly beyond the 
next-nearest-neighbor site, i.e., Cr>2 = for many spin 
models, such as the XY model in transverse field [2], the 
Heisenberg model ground states on ID chains and various 
2D lattices [Q], and even in the Majumdar-Ghosh ( Ji > 0, 
J2 > 0) chain [i^. In the XYZ model in a magnetic 
field, the concurrence extends to several lattice sites, but 
still has finite range [50l |. Long-range concurrences are 
rare; the known examples correspond to specific phase 
transitions Af] or to carefully constructed models and 
geometries [48, 49]. 

In this context, it is noteworthy that our highly frus- 
trated spin rings display nonzero concurrences at all 
length scales, depending on the interaction parameter /3 
(Fig. [2P). For any particular value of (3, only a couple 
of concurrences are nonzero, but not necessarily the ones 
with small r. The r values at which Cr are nonzero, in- 
creases as P changes from to 4, from r = 2 to r = N/2. 
In the region near the ferromagnetic transition /3 = 4, 
there is nonzero concurrence between spins at opposite 
ends of the ring {Cr=N/2 > 0), reflecting the unusual 
long-r ang e quantum correlations of the so-called UDRVB 
state [l9|, [2J] near /3 = 4. The behavior we have found 
indicates that the interaction here acts as a knob to tune 
the distance between spins with nonzero entanglement. 
This tunable distance is quite different from the 'entan- 
glement range' of Ref. [43|; our concurrences are not de- 
caying with distance but are nonzero only at a certain 
distance. 

Fig. [2ji; shows the entanglement between a two-spin 
partition (consisting of site i and site i -\- r) and the rest 
of the ring. Loosely speaking, this entanglement entropy 
is large (small) when the concurrence between the two 
spins is small (large). This reflects the intuitive idea that 
objects highly entangled with each other are generally 
weakly entangled with their environment. 



V. ENTANGLEMENTS AND THEIR 
DERIVATIVES 

Figure [3] displays entanglements 51/2, •S's.i., 5*2/4 for the 
partitions explained in Figure [1] and also the derivatives 
with respect to /3 for the first two cases. 

The block entanglement entropy S1/2 is large com- 
pared to, say, the singlet ground state of the nearest- 
neighbor Heisenberg ground state, in which case 
5i/2 '-1. 28(1.354) for N = 16(20). This is not surpris- 
ing, as the longer-range interaction causes larger entan- 
glement. The Majumdar-Ghosh chain (positive Ji and 
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FIG. 3: 16-spin and 20-spin rings. Top: half-block and sub- 
lattice entanglements, S\/2 (upper dashed), Ss.i. (solid), S2/A 
(lower dashed, magnified), against /? in the singlet region. For 
the N = 20 curves, smaller values of (3 are omitted, to focus 
on the level crossings. Middle and bottom: derivatives of 
Si/2(/?), 5*8.1. (/3), respectively. 



J2) similarly has relatively large block entanglement [2l| . 

The sublattice entanglement 5s. 1. starts from zero at 
/3 = 0, because the two sublattices are decoupled in the 
J2 — > 00 limit. After the first level crossing, 5s. 1. also 
reaches larger values. Note that 5s. 1. is quite large in the 
Heisenberg and Majumdar-Ghosh chains; in comparison 
the ground-state singlets of our system have somewhat 
lower sublattice entanglement. 

The derivatives of 5i/2(/5) and 5s.i.(/3) (Fig. [3] lower 
panels) are similar and have some striking features. First, 
at each level crossing except the leftmost, the deriva- 
tives have shapes reminiscent of a resonance, although 
there is no real divergence. In other words, the entan- 
glement curves S{(3) become steep but not completely 
vertical near each crossing. Second, the S'{(3) curves 
appear to respond (curve upward or downward) already 
some distance away from the level crossing. While phase 
transition language is not completely appropriate for our 
mesoscopic arrangements, the crossings in this sense are 
loosely speaking "second-order" like rather than "first- 
order" like. Finally, the S'{(3) curves at the leftmost 
crossing lacks the resonance-like feature, indicating that 
the first level crossing is different in nature. 

The 52/4 entanglement entropy is peculiarly small 
(shown 4- fold magnified in Fig. [3]), although the rele- 
vant reduced density matrix has the same dimensions as 
the other two partitions. This peculiarity remains unex- 
plained at present. 
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FIG. 4: Fidelity F^{(3) with e = 10"^, and its derivative. Since 
F is generally very close to unity, it is more convenient to plot 
1 — _F on a logarithmic scale. Note that the discontinuity at 
the leftmost level crossing is pronounced in the fidelity but 
almost invisible in its derivative. 



VI. FIDELITIES 

The ground-state fidelity is defined as the overlap be- 
tween ground state wavefunctions at nearby parameter 
values: 



F,(^) = (V'(/3-ie) V(/3+ie) 



Recent studies indicate that the fidelity provides useful 
signatures of quantum phase transitions ^2^]. Since we 
concentrate on mesoscopic spin structures, our purpose 
is not to study phase transitions but to present fidelity 
signatures of the level crossings. 

Figure 0] presents fidelity data for e = 10^'^. The 
derivatives F'{(3) again show striking resonance-like fea- 
tures. The special nature of the leftmost crossing is 
manifested even more dramatically here - the derivatives 
match almost perfectly across this crossing, even though 
the fidelities themselves are discontinuous. 



for adiabatic quantum computing flT'l. It is therefore 
perhaps natural to study the opposite limit of adiabatic 
parameter changes, namely parameter quenches where 
/3 changes as a step function. In an isolated finite spin 
ring, a quench of this type leads to oscillatory behav- 
iors of most c^uantities, including entanglement entropies 
for various partitions; there is no relaxation mechanism 
for the system to reach its ground state. This is shown 
in Fig. [5] (center panel). Studying the evolution after a 
parameter quench from /3i to f3{ involves following the 
wavef unction 

m)) = e^p[~iH{pi)t] 1^(0)) 

explicitly in time, where the initial state \ip{0)) is the 
ground state of H{(3i). We performed this calculation 
by expanding the operator e"'^'^'^'''* to sufficiently high 
order for each time step. 

We also explore relaxation issues, specifically the effect 
of the O-TT-O-TT- crossings, by adding an artificial dissipa- 
tion to the temporal evolution (Fig. rightmost panel). 
Instead of evolving in time t, we evolve in t{l — ij). This 
form of damping ('cooling') preserves the lattice momen- 
tum but not the energy. We display entanglement evo- 
lutions for the large-damping case of 7 = 0.5, and also 
a single example with smaller damping, 7 = 0.1. With 
damping, the sublattice entanglement 5*8.1.(0 relaxes to 
the ground-state 5*8.1. corresponding to the final (3 only 
if the starting state has the same lattice momentum k. 
For N — 12, two level crossings divides the parameter 
space into regions I and III with the same k — n and 
an intermediate region II with fc = 0. The entanglement 
entropy after a quench to /? = 2 (region I) relaxes to the 
'correct' value if the starting ground state is in region I 
or III, but not if it is in region II. 

The structure of level crossings and alternating 
ground-state momenta thus allow selective targeting of 
quantum ground states, and therefore could serve as the 
basis of selective variations of the basic idea behind adi- 
abatic quantum computing. 



VIII. CONCLUSIONS 



VII. TIME EVOLUTION AFTER A QUENCH 

We now turn to the dynamics of entanglement in the 
frustrated ferromagnetic chain. While each entanglement 
measure has its own dynamics when the quantum state 
of the spin ring evolves outside equilibrium, a full study 
of the dynamics of all such measures is outside the scope 
of the present investigation, especially since such com- 
prehensive entanglement dynamics studies are yet to be 
completed for simpler spin systems. We therefore restrict 
the presentation to the dynamics of the sublattice entan- 
glement entropy (5*8.1.). 

The fact that our system has true level crossings, place 
them in the complete opposite limit to what is necessary 



We have identified a mesoscopic system with unusual 
entanglement properties. Our system is the finite-size 
analog of the intriguing ferromagnetic-antiferromagnetic 
chain. The surprising long-distance concurrence we have 
found, with tunable distance, makes this system ideal 
for various quantum information processing tasks. For 
example, entanglement transport is here simply a mat- 
ter of tuning the interaction parameter while remaining 
in the ground state. Our calculations of multi-qubit en- 
tanglement entropies and the fidelities have revealed un- 
expected aspects of the level crossings present in these 
systems. We have also presented entanglement dynamics 
results, again indicating possible quantum information 
processing applications. 
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FIG. 5: Left panel shows sublattice entanglement entropy for 
= 12. Center [rightmost] panel shows evolution of Ss.i. 
after a quench, without [with] dissipation. Initial values are 
A = 0.5 (solid), Pi = 3.5 (dotted), A = 3.85 (dashed). Final 
value is f3{ = 2.0 in each case. Damping constant in rightmost 
panel is shown for 7 = 0.5 for all three A. For A — 0.5, the 
thick solid line has smaller damping, 7 = 0.1. In the presence 
of dissipation, Sa.i. relaxes to the ground state value when 
the initial state has same lattice momentum {k = vr) but not 
when the initial state has different lattice momentum (fc = 0). 



which of the entanglement behaviors reported here also 
carry over to the other systems where such level crossings 
have been observed? To the best of our knowledge, no 
unified study of such crossings exists. Our entanglement 
results call for a general investigation of finite-size level 
crossing sequences occurring as a result of incommensu- 
rate correlations. 

Second, our study of entanglement dynamics treats 
only a minute fraction of the possible kinds of dynam- 
ics, and only a single entanglement measure out of 
many. Since entanglement measures are sensitive to a 
variety of wavefunction characteristics not necessarily 
present in traditional condensed-matter measures such 
as correlation functions, the dynamics of these new mea- 
sures present rich opportunities for characterizing non- 
equilibrium phenomena. Again, comprehensive investi- 
gations are called for, but the study of the evolution 
of various entanglements would be appropriately first 
done in simpler (perhaps exactly solvable) many-body 
systems. 
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